Decreasing Density Gradients in Circumnuclear H II Regions of 
Barred Galaxies NGC 1022, 1326 & 4314 

Jose Franco^, J. Antonio Garcia-Barreto^ and Eduardo de la Fuente^'^ 

Received ; accepted 

version : astro-ph jul 4,2000 



^Instituto de Astronomia, Universidad Nacional Autonoma de Mexico, Apartado Postal 
70-264, Mexico D.F. CP. 04510 



^Facultad de Ciencias, Universidad Nacional Autonoma de Mexico, Mexico D.F. CP. 
04510 



- 2 - 



ABSTRACT 



A re-analysis of the radio continuum emission from circumnuclear regions 
of the barred galaxies NGC 1022, 1326 & 4314 is presented. The spatial 
distributions of Ha and radio continuum are similar and suggest the existence 
of several giant H II regions in the circumnuclear zone. The spectral indices of 
the H II regions, obtained with similar angular resolution and at three different 
wavelengths, indicate that in all three galaxies a2 > —0.1. This result suggests 
the existence of giant H II regions with decreasing density gradients. The 
average density stratifications of these H ii regions could be approximated by 
power-laws of the form Ue oc r~'^, with exponents in the range 1.6 < a; < 2.4. 



Subject headings: galaxies: barred — galaxies: interstellar matter — galaxies: 
giant H ll regions 
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1. Introduction 

The average density profiles of extragalactic H II regions can be derived from 
their thermal radio emission. The thermal radio continuum emission of optically thin 
photoionized plasma with rjj < 1 has a relatively fiat spectrum with Sj, oc z/^°'^, and it 
changes to 8,^ oc when the region becomes optically thick, with r// > 1, and has a sharp 
boundary (e. g. Mezger, Schraml & Terzian 1967; Mezger & Henderson 1967). When the 
electron density varies with distance, say as a power-law with oc r~'^, and the plasma is 
optically thick near the center of the emitting object (e. g. ionized stellar winds or dense 
photoionized molecular clouds), the spectral index varies as S^, oc i/", with —0.1 < a < +2 
for ijj ranging from to oo (Olnon 1975; Panagia & Felli 1975). This behavior is indeed 
observed in stars with extended winds (e. g. Simon et al. 1983) and in some galactic 
ultra-compact and super- ultra-compact H ll regions (Franco, Kurtz & Hofner 2000). In 
particular, the radio continuum emission of the core of galactic source G35.20-1.74 scales as 
Si, oc z/°'^, indicating that its internal density structure is proportional to (Kurtz 2000; 
Franco et al. 2000). 

Here, evidence is given that the average density profile of extragalactic H ll regions 
can indeed be derived from the radio continuum, and discuss the radio emission from the 
circumnuclear regions in three barred galaxies: NGC 1022, 1326 and 4314. Their Ha and 
radio continuum emission at 20, 6, and 2 cm, along with CO, HI, and near IR have been 
reported by Garcia-Barreto et al (1991a,b,c, 1996). In this paper we re-analyze the radio 
continuum data and report the spectral indices of the radio continuum emission from several 
sources in the circumnuclear regions of these galaxies. Their spectral indices between 20, 
6, and 2 cm, with observations at the same angular resolution, indicate contributions from 
both synchrotron and free-free emission. Our results indicate: 1) the existence of several 
circumnuclear H II regions with radio continuum spectral indices > —0.1, and 2) the 
spectral index between 6 and 2 cm from each of the circumnuclear H II regions can be 
explained if the photoionized gas is optically thick and has a decreasing density structure 
with a power-law form Ug oc r~'^, where the exponent is in the range 1.6 < a; < 2.4. 



2. The Ionized Gas in NGC 1022, NGC 1326, and NGC 4314 

NGC 1022 is an SBa(r)p barred spiral galaxy located in the Cetus-Aries group of 
galaxies, at a distance of 18.5 Mpc (Hq = 75 km s~^ Mpc"-^; TuUy 1988). The radio 
continuum emission at 20, 6, and 2 cm (with the same angular resolutions), CO (1-0) 
and CO (2-1) emission, and near IR with a single element detector have been reported by 
Garcia-Barreto et al (1991c). NGC 1326 is an RSBa barred spiral galaxy in the Fornax I 
cluster of galaxies at a distance of 16.9 Mpc (Tully 1988). Radio continuum emission at 20 
cm, 6 cm, and 2 cm, CO (1-0) and CO (2-1) emission, and optical spectroscopy have been 
reported by Garcia-Barreto et al (1991b). NGC 4314 is an SB(a) barred spiral galaxy in 
the Coma I group of galaxies, at an adopted distance of 10 Mpc. Radio continuum emission 
at 20 cm, 6 cm, and 2cm, HI. CO (1-0) and CO (2-1) emission, and near IR with a single 
element detector have been reported by Garcia-Barreto et al (1991a). Images of the optical 
continuum in the I filter (8040 A) and continuum- free Hq;-|-[NII] of these three galaxies have 
been reported by Garcia-Barreto et al (1996). 
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Here we will focus on the properties of the radio continuum emission of the different 
sources in the circumnuclear regions. In all three galaxies the Hq;+[NII] and the radio 
continuum coincide spatially, indicating that the emission is associated with the star 
formation process. The radio continuum emission maps were obtained at the VLA with 
similar high angular resolution beams at 20 cm, 6 cm and 2 cm in order to do spectral 
index analysis of similar regions at the three wavelengths (Garcia-Barreto et al. 1991a,b,c). 
The Ha line emission images were obtained with a CCD camera at the 2.1m telescope in 
San Pedro Martir, B.C. Mexico (Garcia-Barreto et al. 1996). 

Figure 1 is a superposition of the radio continuum emission at 2 cm on the Ho; emission 
of NGC 1022. Both the radio continuum and Ha emitting regions are distributed in a 
similar manner: there is a bright emission region coincident with the compact nucleus, 
and there are two additional compact sources off-nucleus (to the north-east and slightly 
north-west). We have labeled them as NGC 1022:FGF 1 (nucleus), 2 (to the northeast), and 
3 (to the northwest), respectively (see Figure 1 and Table 1). In the original analysis, done 
in 1991, only the two strongest sources, 1 and 2, were considered in the radio continuum 
analysis (Garcia-Barreto et al. 1991c); however, analysis of the Ha distribution suggested 
weak emission from a source to the NW (sec Fig. 1 and Garcia-Barreto et al. 1996), and 
a closer inspection of the radio continuum maps also indicated the existence of source 3 
at the same location. At larger scales. Figure 2 shows that the radio continuum emission 
is found only in the central regions of the disk. The image shows the optical I continuum 
superimposed on the radio continuum at 20 cm. Notice that the position angle of the 
continuum emission between source 1 (the nucleus) and source 2 is almost perpendicular to 
the position angle of the stellar bar {P. A. radio ~ 25° vs. P. A. bar ~ 115). Also, source 3 (NW) 
coincides with a region of distorted isophotes in the optical continuum at P. A. ~ —30°. 

Figures 3 and 4 reproduce the 2 cm radio continuum emission from the circumnuclear 
structures in NGC 1326 and NGC 4314. The source designations are given in Table 
1. The location of these structures indicate that they are most likely the result of gas 
concentrations near (or at) an Inner Lindblad Resonance (ILR) due to the dynamics driven 
by the non-axisymmetric potential of the bar (Schwarz 1984). For instance, for a bar 
angular speed of 36 km s~^ kpc~^ in NGC 4314, the ILR could be at a distance of ~ 450 
pc from the galactic center and, for a bar angular speed of 60 km s~^ kpc~^, the ILR in 
NGC 1326 could be located between 200 and 400 pc from the center (see Garcia-Barreto 
et al. 1991a,b). Gas accumulates near the resonance, massive stars are then formed and 
giant H II regions are created by the strong photoionization field from these newly formed 
massive stars. 



3. Spectral Index of the Radio Continuum Circumnuclear Sources 

We have re-done the maps (see Figs. 1, 3 & 4), and re-derived the spectral indices 
between 20 and 6 cm, and 6 and 2 cm. These indices were computed from the peak fluxes, 
determined from gaussian fits made with the task IMFIT in AIPS. In the case of NGC 1022, 
the linear resolution was 335 by 135 pc and the average spectral index for the compact 
nucleus and NE source between 6 and 2 cm arc a2 ~ +0.12 and +0.15, respectively. In 
contrast, we find af ~ +0.9 for the NW source. For NGC 1326, the linear resolutions are 
295 by 245 pc and the spectral index is ag ~ +0.06 and for NGC 4314, the linear resolutions 
are 195 by 155 pc and the spectral index values are a2 ~ +0.2. The values for all regions 
are given in Table 1. Only the strongest sources, with flux densities at 2 cm larger than five 
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times the rms noise values, are considered. Also, our observations do not resolve individual 
(compact or giant) H ll regions since their average size (less than 100 pc) is smaller than 
our resolution in all three galaxies. 

The spectral index between 20 and 6 cm is negative in all cases, indicating that the 
emission in this wavelength range has a non-thermal origin. In contrast, the spectral index 
from 6 to 2 cm is always positive, suggesting that this emission is dominated by thermal 
bremsstrahlung. This is a consequence of optically thin synchrotron radiation having 
af^ ~ ^0.8 (Niklas, Klein & Wiclcbinski 1997) and thermal bremsstrahlung having a 
between —0.1 and +2, depending on the optical depth (e. g. , Mezger & Henderson 1967). 
Clearly, the values of a2 for all circumnuclear sources in the three galaxies are within the 
range of thermal emission, —0.1 < q;| < +2 (see Table 1). The total radio continuum flux 
at these wavelengths is the sum of synchrotron and free-free emission. However, a detailed 
decomposition of the fluxes obtained with similar beams for NGC 1326 by Garcia-Barreto 
et al. (1991b) indicates a very small synchrotron contribution at 2 cm. A similar exercise 
with the other sources provides the same result. The synchrotron emission probably falls 
more rapidly than the free-free emission rises between 6 and 2 cm, and the positive index 
values are lower limits (but likely close) to the true free-free spectral indices. This is clearly 
seen in Figure 5, that shows vs v for all the sources listed in Table 1: i) the emission 
between 20 and 6 cm falls down less rapidly than typical non-thermal emission (Niklas, 
Klein & Wielebinski 1997), indicating an important thermal contribution at 6 cm, and ii) 
the extrapolation of the non-thermal component at 2 cm is a small fraction of the actual 
emission. 

The lower limits to the average electron density values of the emitting regions are 50 
cm-=^ for NGC 1022 (Garcia-Barreto et al. 1991c), 8 cm-=^ for NGC 1326 (Garcia-Barreto 
et al. 1991b), and 25 cm~^ for NGC 4314 (Garcia-Barreto et al. 1991a). These values are 
similar to the average electron densities, around 10 cm^^, reported for the circumnuclear 
starburst in NGC 1097 (Hummel, van der Hulst & Keel 1987), and to those reported for the 
giant H ll region W3A, within 1 and 80 cm~^ (Kantharia, Anantharamaiah & Goss 1998) 



4. Discussion 

The spectral indices for the circumnuclear regions in our three barred galaxies are 
a\ > —0.1, and they have similar spatial distributions in the Ha, CO and radio continuum 
emission (Garcia-Barreto et al. 1991a, b,c; Benedict, Smith & Kenney 1996). This indicates 
that massive star formation is going on in the giant molecular cloud complexes of the 
observed regions. 

Molecular cloud complexes in our galaxy display highly irregular density and velocity 
distributions, giving the impression of conglomerates of high density condensations, or 
cloud cores, interconnected by a more tenuous intcr-corc medium. These high density 
condensations seem to be the actual sites of massive star formation and they host many 
of the known UCH II and SUCH II regions (e. g. Cesaroni et al. 1999; Kurtz et al. 2000). 
During the formation phase of an H II region, the radiation field creates an ionization 
front that evolves within the density profile of the star forming core. The initial structure 
and early expansion phases of the recently formed H II region, then, are defined by the 
properties of the star-forming core (Franco, Tenorio-Tagle & Bodenhcimcr 1989, 1990). 
Observations of nearby cloud fragments and extinction studies in dark clouds indicate 
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density distributions proportional to r^^, with 00 ranging from 1 to 3 and having an average 
of a; ~ 2 (e. g. Arquilla & Goldsmith 1985; Gregorio Hetem, Sanzovo & Lepine 1988). For 
decreasing density gradients with co > 1.5, the ionization front eventually overtakes the 
shock front and a large region around the core also becomes ionized. All parts of the cloud 
core are set into motion, sometimes driving internal shocks, and instabilities in both the 
ionization and shock fronts generate clumps and finger-like structures (Garcia-Segura & 
Franco 1996; Franco et al. 1998; WiUiams 1999; Freyer, Hensler & Yorke 2000). 

The disruptive effects of photoionization and photodissociation first halt the star 
formation activity (Franco, Shore & Tenorio-Tagle 1994; Diaz-Miller, Franco & Shore 1998) 
and later, once a large fraction of the parental cloud is ionized, accelerate outflows in 
an expanding giant H II region. The combined action of photoionization and the strong 
mechanical energy input from stellar winds and supernova explosions, then, eventually 
destroys the molecular cloud complex and creates large expanding shells (see reviews by 
Yorke 1986, Tenorio-Tagle & Bodcnhcimer 1988, and Bisnovatyi-Kogan & Silich 1995). 
Thus, the non-thermal and thermal radio continuum emission from the circumnuclear 
regions of the three barred galaxies can be viewed as resulting from the stellar energy 
input that is destroying the parental clouds. The negative spectral index between 20 and 
6 cm (except for the compact nucleus of NGC 1022) is likely due to non-thermal emission 
from SN associated to the star-forming activity. The positive index between 6 and 2 cm, 
on the other hand, is likely due to optically thick thermal emission from the expanding 
H II regions. 

The spectra of the free-free emission from an optically thick plasma with various 
decreasing density stratifications have been calculated by Olnon (1975) and Panagia & FeUi 
(1975). The main result of these stTidies is that, for unresolved sources, the spectral index 
is in the range from -0.1 to +2, depending on the particular functional form of the density 
gradient and the optical depth of the region. For a power law of the form Ue oc r~'^, the 
value of the spectral index a, Si, oc u", depends on the exponent cu as (Olnon 1975) 

{2u-3.1)/{u-0.5). (1) 

For a truncated power-law (which removes the pole rig — >■ 00 at r — >^0) and other functional 
forms for the density stratification, the formulae for the spectral index become more 
complicated. For simplicity, and following the results from nearby clouds, the power-law 
case can be used as an approximation to the plasma density stratification. Thus, the 
photoionized circumnuclear clouds detected in the radio continuum can be characterized by 
exponents in the range 1.6 < a; < 2.4 (the values of uj for all sources are listed in Table 1). 
It is important to emphazise, as the referee has pointed out, that these giant extragalactic 
H II regions are usually composed of a collection of individual ionized cores emmbedded in 
a lower brightness extended region. Thus, the density gradients derived from the spectral 
index measurements represent an average density stratification which may be different from 
the actual density profiles of the indiviual ionized cores. 

Summarizing, the radio continuum emission from the circumnuclear regions in the 
three barred galaxies NGC 1022, 1326 and 4314 show a combination of non-thermal and 
thermal features. The fact that all observed sources have similar properties, a negative 
index value between 20 and 6 cm and a positive value between 6 and 2 cm, indicates that 
the radio emission is mainly due to SN and H II regions in massive star forming regions. 
In particular, the free-free emission of the unresolved giant H ll regions indicates density 
gradients shallower than cu ~ 2.5. This result is similar to the one reported by Franco et al. 
(2000) for ultracompact H ll regions, and suggests that many H ll regions, in our Galaxy 
and in other galaxies, are optically thick at radio wavelengths and may provide information 
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of the density stratification of tlie emitting plasma. Tliis was already hinted at the reported 
spectral indices greater than -0.1 in both giant and ultracompact H II regions (Wood & 
Churchwell 1989; Kantharia, Anantharamaiah & Goss 1998), and suggests that further 
analysis of the data may provide valuable information about the properties of star forming 
regions. 
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Table 1: Radio Continuum Spectral Indices 



Galaxy 


Source 


Description 








NGC 1022 


NGC 1022:FGF 1^ 


Compact Nucleus 


-0.83±0.03 


+0.12±0.03 


1.62 




NGC 1022:FGF 2 


North East source 


-0.6±0.1 


+0.15±0.05 


1.63 




NGC 1022:FGF 3 


North West source 


-0.9±0.3 


+0.9±0.3 


2.41 


NGC 1326 


NGC 1326:FGF 1 


Western source 


-0.67±0.03 


+0.06±0.05 


1.58 


NGC 4314 


NGC 4314:FGF 1 


North source 


-0.83lo;o4 


+0.26l°:f 


1.70 




NGC 4314:FGF 2 


North East source 


-0.82±0.2 


+0.23±0.2 


1.68 




NGC 4314:FGF 3 


SoTith East soTirco 


-0.63±0.05 


+0.18±0.1 


1.65 



" a; is the exponent of the density distribution, Ue cx r (see text). 
^ Source also known as NWSS J023832-064039 (Condon et al. 1998) 
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Fig. 1. — Radio continuum emission at 2 cm (contours) superimposed on the Ha emission 
(grey scale) of NGC 1022 (from Garcfa-Barrcto et al. 1991c; 1996). Grey scale: 2.6 to 
170 X 10-^5 ergs cm'^ seg-^/beam. Contours: 2.5, 3, 4, 6, 7, 11, 15, 17, 19, 23, 27, and 31 
times 150 /xJy/beam (rms noise is 150 /xJy/beam). Restoring beam FWHM ~ 3". 7 x 1".5 
at P. A. ~ —11°, corresponding to a linear scale of 335 pc by 135 pc. 

Fig. 2. — Optical continuum emission from the central region of NGC 1022 (contours) 
superimposed on the radio continuum emission at 20 cm (grey scale) (images taken from 
Garcia-Barreto et al. 1991c; 1996). Grey scale: 0.3 to 10 mJy/beam. Contours (in arbitrary 
brigtness units of 10): 9, 9.5, 10, 10.5, 11, 11.5, 12, 12.5, 13, 14, 15, 16, 17, 18, 19, 20, 21, 
22, 23, 24, 25, 27.5, 30, 40, and 50. 

Fig. 3. — Radio continuum emission at 2 cm (contours and grey scale) of NGC 1326 taken 
from Garcia-Barreto et al. (1991b). Grey scale: 0.4 to 2.2 mJy/beam. Contours: -1, 1.5, 3, 
4, 5, 7, and 9 times 250 /xJy/beam (rms noise is 250 /xJy/beam). Restoring beam FWHM 
~ 3". 6 X 3" at P.A. ~ 30° corresponding to a linear scale of 295 pc by 245 pc. 

Fig. 4. — Radio continuum emission at 2 cm (contours and grey scale) of NGC 4314 taken 
from Garcia-Barreto et al. (1991a). Grey scale: 45 to 320 ^Jy/beam. Contours: -2, 2, 2.5, 
3, 4, 5, 5.5, 6, 6.5, and 7 times 45 ^Jy/bcam (rms noise is 45 yuJy/beam). Restoring beam 
FWHM ~ 3". 9 X 3".l at P.A. ~ 90° corresponding to a linear scale of 195 pc by 155 pc. 

Fig. 5. — Sy vs V for all the sources listed in Table 1. For clarity, the sources are grouped by 
galaxy (solid, dotted and dashed lines correspond to sources 1, 2 and 3, respectively), and 
each curve has been multiplied by the factor indicated in parenthesis. 
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